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Abstract In this study, photospheric vector magnetograms obtained with the 
Synoptic Optical Long-term Investigations of the Sun survey (SOLIS), are used 
as boundary conditions to model the three-dimensional nonlinear force- free (NLFF)| 
coronal magnetic fields as a sequence of nonlinear force-free equilibria in spherical 
geometry. We study the coronal magnetic field structure inside active regions 
and its temporal evolution. We compare the magnetic field configuration ob- 
tained from NLFF extrapolation before and after flaring event in active region 
(AR) 11117 and its surroundings observed on 27 October 2010. We compare 
the magnetic field topologies and the magnetic energy densities and study the 
connectivities between AR 11117 and its surroundings. During the investigated 
time period, we estimate the change in free magnetic energy from before to after 
the flare to be 1.74 x 10 32 erg which represents about 13.5% of nonlinear force- 
free magnetic energy before the flare. In this study, we find that electric currents 
from AR 11117 to its surroundings were disrupted after the flare. 
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1. Introduction 

The structure of the Sun's corona is dominated by its magnetic field. To un- 
derstand eruptive events (flares, coronal mass ejections (CMEs) or filament 
eruptions), we need to follow the evolution of the 3D magnetic configuration 



(geometry and topology) (Schrijver and Title, 20111. Knowledge of the amount 



of free magnetic energy and its temporal variation during CMEs/flares will help 



our quantitative understanding of solar explosive phenomena ( |Bleybel et a/., 
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Jing et al, 2010). However, routine measurements of the solar magnetic field 
are mainly carried out in the photosphere. The difficulties of measuring the 
coronal field and its embedded electrical currents thus leads us to use numerical 
modelling to infer the field strength in the higher layers of the solar atmosphere 
from the measured photospheric field. 

Nonlinear force-free field (NLFFF) models are thought to be viable tools 
for investigating the structure, dynamics, and evolution of the coronae of solar 
active regions. NLFFF models were successfully applied to analytical test cases 



(Schrijver et al, 2006 Metcalf et al, 2008), but they were less successful in 



applications to real solar data. Different NLFFF models have been found to 
yield markedly different field line configurations and widely varying estimates 
of the magnetic free energy in the coronal volume, when applied to solar data 



(DeRosa et al, 2009). The main reasons for that problem are (1) the forces 
acting on the plasma within the photosphere, (2) the uncertainties of vector- 
field measurements, particularly of the transverse component, and (3) the large 
domain that needs to be modelled to capture the connections of an active region 
to its surroundings. Therefore, nonlinear force-free modeling is not a routine 
procedure which is guaranteed to produce meaningful results unless the above 



points are taken into account (DeRosa et al, 2009 Schrijver, 2009). In this study, 



we have considered those points explicitly. However, caution is still needed when 
assessing results from this modeling. This is because many aspects of the specific 
approach to modeling used in this work, such as the use of preprocesscd boundary 
data, the missing boundary data, measurement error due to noise, resolution of a 
magnetogram and the departure of the model fields from the observed boundary 
fields may influence the results. 

Solar flares are thought to be powered by the magnetic free energy (i.e., 
the difference between the actual magnetic energy and energy of the equivalent 
potential field) stored in the corona prior to eruption. The storage of free energy 
requires a non-potential magnetic field, and is therefore associated with a shear 



or twist in the coronal field away from the potential, current-free state 1 


Priest 


and Forbes, 2002 Su et al, 2007 


Murray, Bloomfield, and Gallagher, 2011 


). To 



date, the pre-cursors to flaring are still not fully understood, although there is 
evidence that flaring is related to changes in the topology or complexity of an 
active region's magnetic field. 

In this study, we model the coronal magnetic field to determine the sources 
of flaring activity and the temporal evolution of an active region between pre- 
and post-flare stages. Assuming that the evolution of the coronal magnetic field 
above an active region can be described by successive equilibria, we follow in 
time the magnetic changes of the 3D nonlinear force-free (nlff) fields recon- 
structed from two photospheric vector magnetograms taken before and after a 
flare. We use photospheric vector magnetograms as the boundary conditions to 
model the three-dimensional coronal magnetic fields in spherical geometry. This 
enables us to accommodate most of the connectivities within AR 11117 and its 
surroundings. 
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2. Nonlinear force-free field extrapolation 



Except for during eruptions, the magnetic field in the solar corona evolves slowly 
as it responds to changes in the surface field, implying that the electromag- 
netic Lorentz forces in this low-/? environment are relatively weak and that any 
electrical currents that exist must be essentially parallel or antiparallel to the 
magnetic field wherever the field is not negligible. Due to the low value of the 
plasma j3 (the ratio of gas pressure to magnetic pressure), the solar corona is 



magnetically dominated (Gary, 2001). To describe the equilibrium structure of 



the static coronal magnetic field when non-magnetic forces are negligible, the 
force-free assumption is appropriate: 



(V x B) x B = 0, 



(1) 



V - B = 0, 



(2) 



B = B^g on photosphere, 



(3) 



where B is the magnetic field and B ^ g is measured vector field on the photo- 
sphere. Equation ([I]) states that the Lorentz force vanishes (as a consequence of 
J || B, where J is the electric current density) and Equation ^ describes the 
absence of magnetic monopoles. Based on the above assumption, the coronal 
magnetic field is modelled with nonlinear force-free field (NLFF) extrapolation 



(Inhester and Wiegelmann, 2006 Valori, Kliem, and Keppens, 2005| Wiegel 



mann, 2004 |Wheatland, 2004| |Wheatland and Regnier, 20091 |Wheatland and 
Leka, 2011||Amari and Aly, 2010[).~ 



In this study, we extrapolate the three-dimensional NLFF coronal fields from 
the photospheric boundary as successive equilibria. Milne-Eddington inverted 
vector magnetograms, obtained by the the Synoptic Optical Long-term Investiga- 
tions of the Sun survey (SOLIS) / Vector-SpectroMagnetograph (VSM), are used 
as the boundary conditions. From a mathematical point of view, appropriate 
boundary condition for force-free modeling are the vertical magnetic field B n 
and the vertical current J n prescribed only for one polarity of B n (Amari et al., 
|Amari, Boulmezaoud, and Mikic, 1999| |Amari, Boulmezaoud, and Aly, 



1997 



2006). A direct use of these boundary conditions is implemented in Grad- Rubin 
codes (Amari, Boulmezaoud, and Mikic, 1999). Practical computations show, 



however, that the solutions for J n described in DeRosa et al. (2009) B+ and 
B~ can differ significantly for real data containing noise and inconsistencies. 



Wheatland and Regnier (2009) and Wheatland and Leka (20111 implemented 



a scheme which uses both of B+ and B~ solutions together with an error 
approximation to derive a consistent solution. Using the three components of 
B as boundary condition requires consistent magnetograms, as outlined in |Aly| 



( 1989 ) . We use preprocessing and relaxation of the boundary condition to derive 



these consistent data on the boundary. 

For modeling the coronal magnetic field above the active region and its sur- 
rounding, we use the variational principle originally proposed by Wheatland, | 



Sturrock, and Roumeliotis (2000) and later improved by Wiegelmann (2004) 
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in cartesian co-ordinates. The method minimizes a joint measure (C u ) for the 
normalized Lorentz forces (Equation[I]) and the divergence of the field (Equation 
[5]) (each of which should equal zero) throughout the volume of interest, V. 
We have implemented this method for the function (£ w ) in spherical geometry 



(Wiegelmann, 2007 Tadesse, Wiegelmann, and Inhester, 20091 and iterate B to 



minimize . The modification concerns the input bottom boundary field 
which the model field B is not forced to match exactly but we allow deviations 
of the order of the observational errors. The modified variational problem is 



(Wiegelmann and Inhester, 2010 Tadesse et al., 2011): 



B = argmin(£ w ), 



C u — Cc + C j + vC 



f + ^ d + ^ photo , 



(4) 



£ f = / w f (r, d,4>)B- 2 j (V x B) x B| 2 r 2 sin 6drd6d<j), 
Jv 



/ w d (r,0, 


0)|V-B| 


IV 





c 



photo 



(B-B 



obs) 



W(0, 



( B - B obs)^ sin 



where £f and £ d measure how well the force-free Equations (fll) and divergence- 
free ([2]) conditions are fulfilled, respectively. Wf(r, 9, <j>) and w d (r, 6, tf>) are weight- 
ing functions for the force-free and divergence-free terms, respectively, and are 
chosen to be identical for this study. The weighting functions are chosen to 
be constant within the inner physical domain V' and decline to with a co- 



hester 



Tadesse et al. 



sine profile in the buffer boundary region (see Tadesse, Wiegelmann, and In 
(2009, 2011)). The third integral, £photo' * s a sur f &ce 



integral over the photosphere which forces to relax the field B at the photosphere 
towards the measured surface field data, B u g . In this integral, ~W(6,<f>) = 
diag(w rad j a j, w% ans , ^ rans ) is a diagonal matrix which gives different weights 
to the observed surface field components depending on the relative accuracy 
in measurement. In this sense, lacking data is considered most inaccurate and 
is taken account of by setting W(9, <f) to zero in all elements of the matrix. 
SOLIS/VSM provides full-disk vector-magnetograms, but for some individual 
pixels the inversion from line profiles to field values may not have been successful 
and field data there remains undetermined at these pixels. Typically, the field 
is missing where its magnitude is small so that these pixels have a small impact 
on the model even if they were measured correctly. Within the error margin of a 
measured field value, any value is just as good as any other and from this range of 
values we take the value that fits the force-free field best. The different errors for 
the radial and transverse components of are taken into account by different 
values for w rad i a i and u^ rans . In this work we used w ya( \; Yd \ = 20ii^ rariS for the 
surface preprocessed fields as the radial component of B^g is measured with 
higher accuracy. Figure [l] shows the difference in gauss between the measured 
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vector magnetograms and the final values of the model. Hence the final model 
field values on the boundary are consistent with the observed data within noise 
levels (noise due to fluctuations in intensity) which are about 1G and 50G for 
longitudinal and transverse components, respectively. 




1 . . . . 1 . . . 

195 200 


?05 


B/(gouss) 







Figure 1. Surface contour plot of radial magnetic field component and vector field plot of 
transverse field( white arrows) of the difference between original data before preprocessing and 
the final field values after relaxation through the term Cpfo ^ in Equation (14 1 . The maximum 
change in the transverse field is 30G which corresponds to the length of The longest white 
arrow. The vertical and horizontal axes show latitude, 0(in degree) and longitude, </>(in degree) 
on the photosphere. In the area coloured in light blue, field values are lacking. 



We use a spherical grid r, 9, <fi with n r , ng, grid points in the direction 
of radius, latitude, and longitude, respectively. For details of the method, we 



direct the reader to Tadesse et al. (2012). In the present work, we use a larger 



computational domain which accommodates most of the connectivity within AR 
11117 and its surroundings. We also take the uncertainties of measurements in 



vector magnetograms into account as suggested in DeRosa et al. (2009). 



3. Results 

Solar activity on 27 October 2010 was dominated by NOAA AR 11117. A GOES 
CI. 2 flare has been observed at N18°S25° in the active region. NOAA records 
indicate that the event began in soft X-rays (SXRs) which were detected by 
the GOES 15 satellite at 16:59 UT, reaching a peak at 17:04 UT and ending 
at 17:08 UT (see Figure [5]). There were two SOLIS/VSM vector magnetograms 



SOLA: solar_paper_2011.tex; 29 February 2012; 3:06; p. 5 



T. Tadesse et al. 



r ' « ■ ; » » » t ■ ■ * t ■ ' ' i r « 

c-minor flora (class c) 











00:00 04:00 08:00 12:00 16:00 20:00 
Time {UT) 



Figure 2. Solar SXR flux on 27 October 2010 in the wavelength range of 0.1-0.8 nm. Blue 
vertical lines indicate the availability of SOLIS/VSM data. 



taken, one about half an hour before and the other about one hour after the flare. 
In Figure [2j the time the magnetograms were taken is marked by red vertical 
lines at 16:33 UT before the flare and the other at 18:00 UT after the flare. 
As a first step, we remove the net forces and torques from the boundary using 
our spherical preprocessing procedure) Tadesse, Wiegelmann, and Inhester, 2009[ ) 
which brings the photospheric magnetic field closer to the boundary values of 



a force- free field (Molodensky, 1969 Aly, 19891. Then we apply our spherical 



extrapolation scheme using the surface vector field solution obtained from the 
preprocessing scheme. To determine the 3D coronal magnetic field as a nlff 
equilibrium, we need the three components of the magnetic field of AR 11117 
and its surroundings on the photosphere. The computations are performed in a 
wedge-shaped computational box of volume V with 140 x 115 x 250 pixels in 
radial, latitudinal and longitudinal directions. The computational box is large 
enough to include the connectivity between the AR 11117 and its surroundings. 
Figure [3] shows the temporal variations of the photospheric vector magnetic 
components of AR 11117 before and after the flare. The normal component is 
color coded and the horizontal components arc shown as surface vectors. In 
order to quantify the change in the surface vector magnetic field, we computed 
the vector correlation between the fields before and after the flare. We use the 



vector correlation (CVec) (Schrijver et al., 2006) metric which generalizes the 
standard correlation coefficient for scalar functions and is given by 



Cvec — 



(5) 
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(a) Before the flare at 16:33 UT 
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(b) After the flare at 1H:D0 UT 



Figure 3. Surface contour plot of radial magnetic field component and vector field plot of 
transverse field with white arrows. The color coding shows B r on the photosphere. The vertical 
and horizontal axes show latitude, 9(in degree) and longitude, </>(in degree) on the photosphere. 
In the area coloured in light blue, field values are lacking. 
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where Vj and are 2D vectors at grid point i. If the vector fields are identical, 
then CVec = 1; if Vj 1 , then Cyec = 0. The correlation (CVec) of the 2D 
surface magnetic field vectors before and after the flare are 0.96 and 0.87 for the 
radial and transverse components, respectively. From these values we can see that 
there has been some change in the surface magnetic field configuration during the 
flare event. The change in the surface magnetic field is towards an intensification 
in transverse components. This change in transverse components indicates that 
there is also the change in the vertical components of electric current density. 
Figure [4] shows the temporal variation of the vertical electric current density 
on the photosphere. We computed the vertical electric current density on the 
surface using the relation J r = f • V x B (where B the magnetic field) from 
the transverse magnetic field components. We depict a surface color plot of the 
vertical electric current densities before and after the flare in Figures |4^i and 
b, respectively. In order to quantify the change in the radial electric current 
densities on the photosphere, we calculate the total of the absolute value of the 
vertical electric current density before and after flare. We use the pointwise sum 
of the values at nodal points. We find a ratio of the total absolute value of the 
vertical electric current density after the flare to that before the flare to be 0.863. 
The vertical electric current in the active region therefore has decreased after 
the flare. 

Magnetic fields are generally recognized as playing a fundamental role in 
flares. During the course of a flare, the magnetic field is believed to undergo 



major changes (Jing et al., 20091. In this work we have compared the magnetic 
field topologies from two datasets taken before and after the flare using our 
3D NLFF reconstruction method. From a visual inspection of the magnetic field 
lines within the extrapolation volume, we recognize some changes in the magnetic 
field structure during the CI. 2 flare. Figure [5] shows some selected magnetic field 
lines from our reconstructions before and after the flare along with the respective 
SDO/AIA images. Figure [6] shows the magnetic field lines of Figure [5] zoomed 
in and rotated to the limb. In order to compare the fields at the two consecutive 
datasets quantitatively, we computed the vector correlations between the 3D field 
configurations. The correlations (CVec) of the 3D magnetic field vectors before 
and after the flare with respect to the potential field configuration before the 
flare are 0.81 and 0.93 respectively. We see that the magnetic field configuration 
after the flare looses some of its non-potentiality. 

The energy stored in the magnetic field as a result of fieldline stressing into a 
non-potential configuration has been identified as the source of flare energy. The 



study by Jing et al. (2010) confirms that there is physical link between magnetic 
energy and flare occurrence in active regions. Study of the temporal evolution of 
the free magnetic energy indicates that it varies before and after the flare events 



(Jing et al, 2009). There is strong need to estimate this free energy numerically. 
One way to estimate the energy budget of active regions is to reconstruct the 
three-dimensional (3D) coronal field from the measured photospheric boundary 
based on the force-free assumption. We compute the free magnetic energy from 
the excess energy of the extrapolated field beyond that of the potential field 
which satisfies the same B n u„-r boundary condition. Similar estimates have been 



labal 

made by Regnier and Priest|(|2007[),|Thalmann, Wiegelmann, and Raouafi|(|2008[), 
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Figure 4. Surface contour plot of the radial component of electric current density. The color 
coding shows J r on the photosphere. The vertical and horizontal axes show latitude, 6 (in 
degree) and longitude, <j> (in degree) on the photosphere. In the area coloured in light blue, 
field values are lacking. 
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(a) AIA image at 16:33 UT 





(b) NLFF solution at 16:33 UT 




(c) AIA image at 18:05 UT 



(d) NLFF solution at 18:00 UT 



Figure 5. SDO/AIA ( textitSolar Dynamics Observatory/ Atmospheric Imaging Assembly) 
17lA images and their respective selected magnetic field lines reconstructed from SOLIS 
magnetograms using nonlinear force-free modelling. 



and Tadesse et al. (2012) for active regions observed at other times. From the 
corresponding potential and force-free magnetic field, Bp t and B, respectively, 
we can estimate an upper limit to the free magnetic energy associated with 
coronal currents 



EfTve. — E. 



nlff " ^pot 



E-^ni — 



1 

8tt 



/ { B l\S ~ B pot) r 2 sin9drddd(f>. 



(6) 



The computed energy values are listed in Table [T] The change in the free energy 
during the flare is about 1.74 x 10 32 erg. The magnetic energy of the potential 
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(a) NLFF solutiuti before the flare at 16:33 UT 




(b) NLFF solution after the flare at 18:00 UT 



Figure 6. Some magnetic field lines reconstructed from SOLIS magnetograms. These Figures 
are obtained by zooming and rotating Figures [5J5 and d to the solar limb. 
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Table 1. The magnetic energy obtained from potential and NLFF field extrapola- 
tions before and after the flare within the computational box. 



Events 


E nlff (10 32 erg) 


£ po t(l° 32 erg) 


£ frcc (10 32 erg) 


Before the flare at 16:33UT 


12.89 


10.84 


2.05 


After the flare at 18:00UT 


10.93 


10.62 


0.31 



field configuration is about 10 x 10 32 erg. i? m fj exceeds £p t by only 15.9% and 
2.8% for before and after the flare, respectively. To estimate the uncertainty in 
the numerical result, the code was applied to the original SOLIS data to which 
artificial random noise had been added in the form of a normal distribution 
with an amplitude of approximately 1 G in the longitudinal and 50 G in the 
transverse component. The chosen noise amplitudes are based on the sensitivity 
of the VSM instrument. Hence, we found that the evaluated relative error of the 
energy estimation is about 0.4% for the potential and 1% for the NLFF field (i.e. 
E pot ±0.044 x 10 32 erg and £ nlff ±0.129 x 10 32 erg, respectively). The available 
free magnetic energy is approximately 10 32 erg with a relative error of about 4% 
(i.e. £'f roo ±0.082x 10 32 erg). In principle high spatial resolution of a magnetogram 
gives a higher magnetic energy because small-scale magnetic variations are better 
resolved. How much better resolution contributes to the total energy estimate 
depends on the steepness of the spectral energy distribution. Steeper spectra 
yield less additional energy if the resolution is enhanced, flatter spectra yield 
more additional energy. In addition to our relative error estimates above, low 
spatial resolution might increase the error. 

To study the influence of the use of preprocessed boundary data along with the 
departure of the model field from the observed boundary fields on the estimation 
of free-magnetic energy, we have computed the magnetic energy of the potential 
field and the NLFFF obtained from the original data without preprocessing and 
with preprocessing. As the preprocessing procedure filters out small-scale surface 
field fluctuations, the magnetic energy of NLFFF obtained from preprocessed 
boundary data is smaller than the corresponding energy without preprocessing. 
The energy of the potential field obtained from boundary data with and without 
preprocessing are close in value, since the potential field calculation makes only 
the use of the radial magnetic field component which is not affected too much 
by the preprocessing. The magnetic energy computed from the original data 
without preprocessing is about 13.67 x 10 32 erg which is about 6% higher than 
the one obtained from preprocessed and modified observational boundary data. 
However, this energy does not correspond to the nonlinear force-free magnetic 
field solution since the original boundary data without preprocessing is not a 
consistent boundary condition for NLFFF modeling. 

In our previous work (Tadesse et al., 2012), we have studied the connectivity 
between three neighbouring active regions. In that study, we investigated the 
three ARs were found to share a significant amount of magnetic flux compared 
to their internal flux connecting one polarity to the other. In terms of the electric 
current they were much more isolated. In this work, we study the connectivity 
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Blue = To the surrounding 

Red= Within AH 

Green= Reach top boundary 

(open) 




Figure 7. The connectivity between AR 11117 and its surroundings. The black rectangle 
shows the domain of AR 11117. The color coding shows B r on the photosphere. The red field- 
line represents fieldlines connecting opposite polarities within the active region. The blue one 
represents fieldlines connecting the active region with its surroundings both on the photosphere 
and the side boundaries (categorized as elsewhere). The green one shows fieldlines, which leave 
the computational box though its top boundary and may be considered to be "open". 

between AR 11117 and its surroundings before and after flare. In order to quan- 
tify these connectivities, we have calculated the magnetic flux and the electric 
currents shared between the active region and its surroundings. For the magnetic 
flux, e.g., we use 



where the summation is over all pixels of AR a from which the field line ends in 
AR^j or i £ AR Q || conjugate footpoint(i) e ARp. The indices a and (3 may take 
values between 1 and 3. The index number 1 corresponds to AR 11117, index 
number 2 to its surroundings and the side boundaries and index number 3 to the 
top boundary (see Figure [7]). For the electric current we replace the magnetic 
field, Bi • f, by the vertical current density Jj • r in Equation Both Table [2] 
and [3] show the percentage of the total magnetic flux and electric current shared 
between the AR 11117 and its surroundings before and after the flaring event. 
For example, first column of Table [2] shows that 82.23% of positive polarity of AR 
11117 is connected to negative polarity within itself; line 2 shows that 41.51% of 
positive/negative polarity of AR 11117 is connected to positive/negative polarity 
of its surroundings including the side boundaries of the computational box, and 
line 3 shows that 74.78% of the total magnetic flux of top boundary of the 
computational box connected to the positive/negative polarity of AR 11117. 
Table [3] shows the electric current connectivity we have calculated applying the 
same technique. Figure [4] shows that the vertical electric current density has 




(7) 
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Table 2. The percentage of the total magnetic flux shared between AR 11117 
and its surroundings. <J>ii, $22 and $33 denote magnetic flux of AR 11117, 
outside AR 11117 on the photosphere including the side boundaries and the 
top boundary of the computational box (see Figure n\, respectively. 



Before the flare After the flare 



<5> a p 


= 1 


2 


3 


= 1 


2 


3 


a = l 


82.23 


12.61 


5.16 


88.37 


7.36 


4.27 


2 


41.51 


50.94 


7.55 


35.18 


56.24 


8.58 


3 


74.78 


25.22 


0.00 


65.49 


34.51 


0.00 



Table 3. The percentage of the total electric current shared between AR 
11117 and its surroundings. In, i22, and ^33 denote electric current within 
AR 11117, outside AR 11117 on the photosphere including the side boundaries 
and the top boundary of the computational box (see Figure [7|, respectively. 

Before the flare After the flare 



la/3 

a = l 
2 
3 



13 = 1 


2 


3 


13 = 1 


2 


3 


93.02 


6.98 


0.00 


98.06 


1.94 


0.00 


34.41 


65.59 


0.00 


26.38 


73.62 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 



decreased after the flare. In this study (see Table[3]), we found that AR 11117 is 
even more isolated in electric current from its surroundings after the flare. It is 
noteworthy that modeling this active region with small cartesian box would lead 
to wrong NLFF model solution as there are currents crossing its boundaries. 

4. Conclusions 

We have investigated the coronal magnetic field associated with AR 11117 and its 
surroundings observed on 27 October 2010 by analysing SOLIS/VSM data. Two 
vector magnetograms with a time cadence of 1 hour and 28 min were available to 
investigate the magnetic energy content of the coronal field during the CI. 2 flare 
observed by GOES. We have used an optimization method for the reconstruction 
of nonlinear force- free coronal magnetic fields in spherical geometry by restricting 



the code to limited parts of the Sun ( 


Wiegelmann, 2007 


Tadesse, Wiegelmann, 


and Inhester, 2009 


Tadesse et al, 2011 


Tadesse et al, 2012 


1. 



We have studied the time evolution of the magnetic field from before to 
after the flare. We found that there is some rearrangement in the magnetic 
field configuration after the eruption. The magnetic energies calculated in a 
large wedge-shaped computational box above the active region and its sur- 
roundings decreased after the flare, indicating that the field looses some of its 
non-potentiality. However, caution is needed when estimating the free magnetic 
energy using NLFFF modeling. Many aspects of the specific approach used in 
NLFFF modeling may influence the results. This is the first study which contains 
AR and its surroundings with a flaring event in our model. It was made possible 
by the use of spherical coordinates and allowed us to analyse connectivity be- 
tween AR and its surroundings. Modeling an active region in cartesian geometry 
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would lead to wrong NLFF model solution as there are currents crossing the 
small cartesian box enclosing it by excluding its surroundings. We propose to sys- 
tematically study the effect of using cartesian box over spherical wedge-shaped 
box for NLFF solutions. 

High cadence magnctogram observations are necessary when we study the 
magnetic field topology and energy variations associated with CME/flare erup- 
tions. In this sense, it is worth mentioning that the Helioseismic and Magnetic 
Imager (HMI) on board the Solar Dynamic Observatory (SDO) is the first 
instrument to provide routine measurements of the full-disk photospheric vector 
magnetogram data with high spatial and temporal resolution under seeing-free 
condition. We anticipate extending the current study with soon to be released 
SDO/HMI full-disk vector magnetograms. 
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